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ABSTRACT 

Angular momentum dependence of nuclear level densities at finite tem- 
peratures are investigated in the static path approximation(SPA) to the par- 
tition function using a cranked quadrupole interaction Hamiltonian in the 
following three schemes: (i) cranking about x-axis, (ii) cranking about z-axis 
and (iii) cranking about z-axis but correcting for the orientation fluctuation 
of the axis. Performing numerical computations for an sd and a pf shell 
nucleus, we find that the x-axis cranking results are satisfactory for reason- 
ably heavy nuclei and this offers a computationally faster method to include 
the angular momentum dependence at high temperatures in the SPA ap- 
proach. It also appears that at high spins inclusion of orientation fluctuation 
correction would be important. 
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1. Introduction 



The level density is one of the important physical quantities appearing in 
the statistical analysis of nuclear reactions It also provides a reasonable 
basis for testing the applicability of the various approximations used for many 
particle systems. The most commonly used mean field approximations (MFA), 
like Hartree-Fock|2[], Hartree-Fock Bogoliubov|3|, M have been successful in 
describing the structure of the nuclei upto a very high spin state along the 
yrast line as well as for the low-lying excited states. However, as one moves 
far away from the yrast line, the theoretical predictions based on the MFA 
show an inconsistency with the experimental data on the shape transitions]^], 
strength function for giant dipole resonance(GDR)||, angular distribution 
of GDR 7-rays0 etc. A simplest way to improve upon this inconsistency 
is to account for the statistical fluctuations associated with various degrees 
of freedom (e.g. pairing and qadrupole degrees of freedom) quadrupole and 
pairing degrees of freedom) caused by the symmetry breaking [g] required to 
incorporate the correlations. As a first step one includes the effects of fluctu- 
ations by considering an isothermal probability distribution as a function of 
an appropriate set of phase variables. In this way one finds a remarkable im- 
provement in understanding the experimental data. However, Alhassid and 
Bush have shown recently |J that this method of including statistical fluctu- 
ations in the MFA leads to a serious problem of overestimating the values of 
the level density. On the other hand, they find that the static path approx- 
imation(SPA) |K| H, [12] which provides a natural framework to deal with 



the statistical fluctuations yields an almost exact value of the level density 
at high temperatures. 

The SPA has been applied for the model studies of the nuclei at finite 
temperatures. All these model studies indicate that SPA becomes almost ex- 
act at high temperatures. Furthermore, the results at low temperatures can 
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be improved by including the contribution of Fock term[|T3] or the quantum 
correlation corrections in the random phase approximation(RPA) method 
14, [15|]. For a realistic application, however, one must restore the broken 



symmetries in order to characterize the system by a set of physically observ- 
able quantum numbers, e.g., angular momentum, particle numbers etc. In 
ref. [IB], the usual exact projection method has been employed to conserve 
the particle number at finite temperature. Recently, this method has been 
extended to project out a general symmetry from the partition function in- 
volving a one-body statistical operator [[17]]. However, this method is rather 
involved to be numerically tractable for heavy nuclei. For instance, the pro- 
jection of an exact angular momentum (i.e. SO (3) symmetry group) state 
would require an additional three dimensional integration besides those over 
the phase variables appearing in the SPA representation of partition func- 
tion^. g. see eq. (0)). On the other hand, by applying SPA to a standard 
J^-cranking Hamiltonian, in refs. [ TBI [0| we have studied the thermal prop- 
erties of a medium heavy rotating nucleus, 64 Zn. It is shown that the results 
for the energy, level density, moment of inertia etc. provide a reasonable 
physical insight at high spins and temperatures. In the present work we 
further investigate the angular momentum dependence of the level densities 
at finite temperatures in SPA using cranked quadrupole Hamiltonian in the 
following three schemes: (i) cranking about x-axis as just mentioned, (ii) 
cranking about z-axis and (iii) cranking about z-axis but correcting for the 
orientation fluctuation of the axis. 

In the next section we present the basic expressions for computing the 
level densities in different cranking schemes mentioned above. Section 3 
contains the numerical details and discussions on the results obtained for an 
sd ( 24 Mg ) and a pf ( 6A Zn ) shell nucleus. Finally, a brief summary and 
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conclusions will be presented in section 4. 



2. Theoretical framework 

The grand canonical partition function for a cranked Hamiltonian is given 

as: 

Z(fi p , fi n , U, T) = rre -(^-MP^.-M»JVn- W Ji)/T (1) 

where i denotes the cranking axis in laboratory frame, the trace is taken over 
all the many-body states of the system described by the Hamiltonian H. 
The quantities fi p , \i n and ui are essentially the Lagrange multipliers used to 
constrain the average number of protons N p , neutrons N n and the component 
of total angular momentum along the cranking axis, i.e, Jj. The Hamiltonian 
H considered here is 

H = H -\xY,Q»Ql ( 2 ) 

where, H Q = Ylih (i) stands for the spherical part of the Hamiltonian with 
h (i) corresponding to the basis state single-particle (sp) energies, x is the 
quadrupole interaction strength and = r 2 !^ is the quadrupole moment 
operator. 

In SPA, the eq.(|I|) can be written as follows, 

Z(fjL p , [i n , u, T) = TrV(/j, p , fi n , u, T) (3) 
where, V is a one-body static path statistical operator as given below, 

V = J V{a) e -(H'(v)-VpNp-tinN n -u,jz)/T ^ 

where, / T>(a) denotes the integration over the temperature independent ( or 
static) paths in the space of auxilliary fields, a = {ai, cr 2 , ...} and H'(cr) is a 
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one-body Hamiltonian. Using ref. pi we have 



/ V{a) = ( 2^T )5/2 / n^-*^' 2 (5) 

^^.-xEH)"^ ( 6 ) 

The above eqs. (|5]) and (Q) can be represented conveniently in the intrinsic 
coordinate system using the eqs. and (§) given below, 



1 '■'^n„I);, l ,(o.0.c) (7) 



X 



where, a± 1 = 0, a± 2 = ftsiwy/yl and a = Pcosj and D is the standard 
Wigner's function. 

Ji^KMiW) (8) 
j 

where, R is the transformation matrix which essentially comprises of usual 
rotation D matrices and the indices % and j denote the component in the lab 
and the intrinsic frame, respectively. Let us first consider the third scheme 
(as discussed above). Using eqs. (|3|) - (§), we have (see ref. flOU ) 

cy r ir. ft™ r n rPmax 
Z = 2ir{-—f /2 dip dOsinO d(3 /3 4 
2irT Jo Jo Jo 

x f /3 d 1 | sm3 7 | e -^Tre-^'^-^ )/T (9) 

J o 

In the above a = M u) ° . M is the nucleon mass,x = 70A~ 1A MeV taken from 

X 

ref. |2]] , cj = 41/A 1 / 3 MeV and # /£J = Ei 



= /i - Mu; 2 r 2 /?[c s 7 Y 2 ,o + ^m 7 (^2,2 + Y 2 ,- 2 )] 

— djjcos^jx — sin6cosipj x + sinOsinipjy} (10) 
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The partition function for other two cases can be obtained simply by replacing 
the coefficient of u (curly bracket) in eq. fllOf) by j x or j z . We now give the 
basic expressions used in computing the energy and level density in terms of 
the partition function for case (ii) and (iii). 

E = T 2 ^- + p p N p + p n N n + ooM (11) 

where, 

^dlnZ , , 

N '^ T a^ (12) 

and 

^dlnZ , , 

M = T— (13) 

is the component of the total angular momentum along the z-axis. 

The level density p(E, J) at a fixed energy and angular momentum is 
given byf22|, [23| 



p(E, J) = p(E, M — J) — p(E, M = J + 1) « - M g M M) \m=j + i (14) 

where, p(E, M) is an inverse Laplace transform of the partition function 
which can be written in a saddle point approximation as 

»(B,M) = (15) 

D is the determinant of a 4 x 4 matrix with elements 
d 2 lnZ 

dij = ^ q x ! x i = (V T > ~lh/ T > -Vn/T, -oj/T) (16) 

On similar lines, the level density for the case (i) can be obtained directly 
using eq. (|I5D with the partition function and its derivatives evaluated at a 
fixed value of J (see ref. [[IIJ for detail). 

Finally, we give expressions for the level density parameter a. Most com- 
monly one uses, 

E* = aT 2 (17) 
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where, E* = E(T)—E(T = 0) is the excitation energy at a given T. However, 
we know that SPA is not applicable in T — > limit. We therefore use, 

S = 2aT (18) 

where, S — (E — F)/T is the entropy with F being the free energy. 
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3. Numerical details and results 



3.1 NUMERICAL DETAILS 

In this subsection we shall present a detailed description ot the numerical 
procedure for the calculation of level densities in different cranking schemes 



as outlined in the previous section. Some detailes as presented in ref. |L9 
are not reported here. 

We see from the eqs. (|TTJ) , QT2D, (|TJP and ( P~5| ) that the evaluation of the 
physical quantities, like energy, particle numbers and level density require 
the computation of partition function (|9]) and its first and second derivatives 
with respect to T, fi PjTl and uj. To begin with, in a suitable model space (e.g. 
see table 1), we diagonalize the one-body Hamiltonian ( |T(]| ) at mesh points 
in the space of quadrupole degrees of freedom, i.e, /3, 7, 6 and ip (note that 
eq. (U) is independent of the variable 0). Now, the trace appearing in eq. @ 
can be obtained for a given value of [i p ^ n and uj using the following equation 



(see also [[□], [LJ 



where, e-s are the eigenvalues of the Hamiltonian (|TUD and the index i runs 



over all the deformed single particle(sp) orbits for proton as well as neutron. 
Equation (§) is then computed by performing a numerical integration using 
eight Gaussian points for each of the variables and taking f3 max = 0.5. Having 
determined Z as a function of [i p , n and uj at a fixed T, we calculate its first 
derivatives with respect to \i P)n and uj and adjust these such that the eqs. 



(|12D and (|13|) are satisfied for the desired values of N Pjn and M, respectively. 
The energy can be easily calculated using eqs. (|TT|) . 

Now, we come to the numerical evaluation of the level densities for three 
different cases discussed in section 2. For the case (i), i.e J^-cranking, we 



8 



first compute the partition function for a desired value of J by adjusting u 



such that the r.h.s. of eq. fll3|) leads to < J x >= J J(J + 1) for a given J. 

The level density p(E, J) is then evaluated using this partition function in 
eq. (|15|), where E = E(T) represents the average energy at a temperature 
T. Note here that for the case (i), the integration over 9 and ip appearing in 
the eq. (^) can be replaced by a factor Att. For the case (ii), i.e J^-cranking, 
it is apparent from the eq. ( |14|) that the quantity p(E, J) can be calculated 
simply by taking the difference between p(E, M = J) and p(E, M — J + 1) 
at a fixed value of E. This means that the temperature has to be adjusted 
such that the energy of the system remains the same at M = J and J+l. As 
an illustration, we show in figure ffl, the variation of energy with temperature 
for M = 2, 3, 4 and 5 for 24 Mg. The points A and B on the horizontal 
dashed line indicate the change in the temperature (~ 0.3 MeV) with M 
changing from 2 to 3 at E = -47.5 MeV. It may be noticed from the figure 
that the temperature difference for a given E in the high temperature region 
decreases. However, instead of adjusting T, we adopt a slightly different 
but relatively faster method which is as follows. To illustrate the method, 
let us take the calculation of p(E, J = 2). We first calculate the values of 
p(E(T),M) for M = 2 and 3 at the temperatures 0.4, 0.6, 0.8,. ..,4.0 MeV. 
Then for each value of p(E(T),M = 3) we calculate p(E(T'), M = 2) using 
standard interpolation method|24j] such that E M=3 (T) = E M=2 (T"). A similar 
numerical procedure is followed for the evaluation of level density in case 
(ii). We must add here that the case (iii) requires the diagonalization of a 
complex one-body Hamiltonian, since J* = —J y , unlike the other two cases. 
However, the problem of diagonalizing anxn complex Hermitian matrix can 
easily be mapped into a diagonalization of a 2n x 2n real symmetric matrix. 




0- 
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3.2 RESULTS AND DISCUSSIONS 



Following the procedure as described above, we have performed the nu- 
merical calculations for an sd ( 2A Mg) as well as a pf shell ( 64 Zn) nucleus. 
In this subsection we shall present the results for level densities and inves- 
tigate its angular momentum dependence at finite temperatures using SPA 
to three different cranking schemes as discussed in the earlier sections. We 
must mention that the model spaces (see table 1) considered here are pre- 
cisely the same as used in refs. |25j and |TTJ for the study of 24 Mg and 64 Zn 
at finite temperatures, respectively. From table 1 it is clear that for 2A Mg 
there are 8 valence particles with 16 as a core and for 6A Zn there are 24 
valence particles with 40 Ca as a core, 
(i) ™Mg 



The level density 



Usually the level density is estimated through the level density parameter 
a with numerical value between A/8 - A/10 given by the Fermi-gas model. 



The experimental data|26|, ^7], ^8| suggest that for A ~ 160 the parameter a 
decreases from A/8 at low temperatures to A/13 at T ~ 5 MeV. Using J x - 
cranking we have studied here the temperature as well as angular momentum 
dependence of the inverse level density parameter K{= A /a). Figure |2] shows 
the variation of K with temperature for a few values of angular momenta, 
J = 0, 2, 4 and 6. We find that the level density parameter a decreases 
with temperature and angular momentum. On the average, the value of a 
decreases from A/6 at T = 1 Mev to A/12 at T = 3.0 MeV. More on this 
will be discussed below. 

In figure ^ we have displayed the J x cranking results for the level density 
as a function of energy for a few values of angular momenta as indicated. 
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Similar plots are shown in figures (^) and (|5]) where the angular momentum 
dependence is brought through cranking along the z-axis without and with 
orientation fluctuation corrections(OFC), respectively. We find from these 
figures that the qualitative behaviour for the level densities obtained in dif- 
ferent cranking schemes considered here is more or less the same. To fascilate 
further discussion we denote the level densites obtained using (i) J^-cranking, 

(ii) J^-cranking and (iii) J z cranking with OFC by p x , p z and p z , respectively. 
We now analyze the level densities for J = 0. It is well known that in J x - 
cranking, J = simply means uo = (or no cranking). So, strictly speaking, 
the quantity p x (E, J = 0) essentially corresponds to the intrinsic or unpro- 
jected level density. On the other hand, using eq. flHD, for the case (ii) and 

(iii) we have, p(E, J = 0) = p{E, M = 0) - p(E, M = 1). This implies that 
for J = 0, Px > p 3tZ (since p z (E, M = 0) = p z (E, M = 0) = p x (E, J = 0). 
In addition to this we also note that p z (E, J = 0) ~ p z (E, J = 0) and the 
difference, p xz = p x (E, J = 0) — p z (E, J = 0) is quite large. Moving towards 
higher values of J, we find that p z < p z for J = 4 and 6 and p x > p z for 
all J. However, the difference p x — p z reduces with increase in J. For exam- 
ple, we find that p xz = 2.718 , 2.459, 2.225 and 1.234 for J = 0, 2, 4 and 
6, respectively, at E = -40.0 MeV. For J = 6, the value of p xz is obtained 
after smoothing out the fluctuations in Inp by an exponential fitting. An 
interesting feature we have noticed is that for a fixed value of E and J the 
level densities p x , p z and p z do not correspond to the same temperature. As 
we pointed out in section 3.1, for the computation of p z and p z one has to 
adjust the temperature in such a way that 

p(E(T), J) = piEin), M = J) - p(E(T 2 ), M = J + 1) (20) 

where, T\ and T 2 satisfies, 

Em=j(T x ) = E M=J+1 (T 2 ) (21) 
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We should note that T appearing on the left hand side in eq. |20] is not a 
quite well defined quantity. Just for the comparison if we define an average 
temperature T av = (7\ + Tz)/2, we find that the T av in p z and p z for a 
fixed E and J are quite close. However, looking into the figure |l| we can 
say that the T av is quite different from the corresponding T in p x (E, J) at 
low temperatures. Let us now turn towards the discussion on the unique 
feature for the level densities corresponding to J = 6. Figures § and f| show 
large fluctuations in the values of p x and p z , respectively, for J=6. We find 
that it is quite difficult to adjust the value of u in the r.h.s. of eq. 



for < J x >= y6(6 + 1) as well as M = 7 upto a desired precission of 10 4 . 
In principle, one may achieve this precission by varying u in a very small 



steps. However, it must be clear from eqs. (0) and (fLOP that for computing 
the partition function, each time for a fixed value of u and T one has to 
diagonalize the one-body Hamiltonian h luj at mesh points in the space of 
quadrupole variables (see sec. 3.1). So, with a smaller increment in u it 
would take enormous computer time to reach a desired value of < J x > and 
M. We therefore fix the accuracy to be ~ 10 _1 . The fluctuations in the 
values of the level densities p x and p z are indicative of numerical difficulty 
in adjusting u for J = 6. For example, the points a and b in the figure |3] 
correspond to the angular momentum J a = 6.15 and J& = 5.85, respectively, 
(ii) 64 Zn 



The level density 

Figure |6| shows the dependence of inverse level density parameter K on 
the temperature and angular momentum. The curves in this figure essentially 
depict a similar trend as that for the case of 2A Mg in Fig. 2, i.e, a decreases 
with temperature and angular momentum. However, one may note that the 
values of K are higher for 64 Zn. This feature is consistent with a recent 
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theoretical ||29|| investigation on the temperature and mass dependence of the 
level density parameter in a realistic model. Furthermore, we observe from 
the figures ^| and |6] that K increases almost linearly with temperature even 
at low T. On the other hand, in ref. [29] it is shown that for A = 40 and 
60 it remains more or less a constant up to T~ 1.5 MeV and then increases 
slowly. It is clear from the ref.|50| that K becomes almost a constant at 
low temperatures only when the frequency dependent effective mass (m^) is 
introduced to simulate the effects of collectiveness. At high temperatures, 
the influence of effective mass disappears, i.e m^/m approaches unity. This 
reveals the fact that at low temperatures, the pairing correlations and RPA 
corrections should be included in the present approach. It must be mentioned 



here that a similar plot for the parameter K has been shown in ref. JJTJ for 
im Er using a finite temperature mean field approximation. It is found that 
the value of K is independent of temperature as well as angular momentum 
for 0.5 < T < 1.4 MeV. At low temperature, T < 0.5 MeV there is a 
decrease in K at all spins. However, the variation of K for T > 1.4 MeV is 
quite close to as shown presently in figures || and §. 

In figures [7| and |8| we have displayed the level densities p x and p z varying 
with energy at fixed values of angular momenta. As in the case of 2A Mg, here 
too we find that the J^-cranking yields higher values of the level density, i.e. 
p x > p z . For J = 0, difference between the values of p x and p z is very 
small compared to that for J / 0. However, the percentage difference, i.e 
(Px — Pz)/ Px is smaller for 64 Zn. 

So far we have only studied the behaviour of the level denstities obtained 
by applying the SPA in different cranking schemes. Apart from the fact 
that SPA is not applicable at very low temperatures, the present cranking 
approach shows limitations at low spins also. For this purpose we com- 
pare the Jj.-cranking results for level density in 2A Mg with the one reported 
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recently [T7|] for 20 Ne using an exact angular momentum projection within 
SPA. As expected, we find the large difference between the results for J = 
0. However, for J > 4, the present results seem to be reasonable. On the 
other hand one would expect that the level density based on other cranking 
along z-axis should be quite realistic even at low spins. It is surprising that 
the value of the level densities for J = are the highest even in the schemes 
(ii) and (iii). This may be the indicative of the necessity of an exact angular 
momentum or K projection at low spins. Finally we would like to add that 
the level density versus energy plots discussed above provide only an easy 
means to compare the results obtained using different schemes. However, for 
a more realistic comparison one needs to study the variation of level density 
with the excitation energy. For example, in figure |^ we have plotted p x vs 
excitation energy E* (obtained from the use of eqs. ([18]) and ([17])) for fixed 
values of angular momenta. Unlike in figures [7| and |8], we find that for E* = 5 
- 10 MeV (about netron separation energy i.e. very low T), the level density 
at J = 16 is slightly higher than that for J = 8. But at higher energies (T>1 
MeV) J = 16 curve falls below that for J = 8 and there is a tendency of 
saturation for all the curves. 

4. Summary and conclusions 

We have investigated the angular momentum dependence of the level density 
of 24 Mg and 64 Zn using static path approximation method with a quadrupole 
interaction Hamiltonian in three different cranking schemes: (i) cranking 
about x-axis, (ii) cranking about z-axis and (iii) cranking about z-axis in- 
cluding orientation fluctuation corrections (only for 24 Mg). In a comparative 
study we make the following observations. 

For 24 Mg we find that qualitatively p z = p z for J = and 2, with p z < 
p z for J = 4 and 6. This may be an indication of the importance of the 
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orientation fluctuation at high spins. However, we also notice that at high 
energy, e.g. E(T) > -42 MeV even at J = 6 p z ~ p z . On the other hand 
Px > Pz for J < 6 and p x ~ p z at J = 6. As expected, this shows that x-axis 
cranking may be a good approximation at high spins only. 

Now coming to the relatively heavier nucleus 64 Zn we find that p x « p z at 
all spins J = - 28. Thus for a heavy nucleus and particularly at high spins 
the x-axis cranking may provide a good prescription for the computation of 
spin dependent level densities. Due to computation limitatins we have not 
yet included the orientation fluctuation correction for 64 Zn. We propose to 
investigate on this in our next calculation. We prefer to incorporate orienta- 
tion fluctuation corrections in the x-axis cranking scheme as it also provides 
us with an opportunity to compute moment of inertia as S =< J x > /uj. 



Recently [19| we hve studied the variation of function of spin and 

temperature for 64 Zn. 
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Table 1: Single particle energies (in MeV) for the spherical orbits present in 
the model space for 24 Mg and &A Zn. 



Model space for 2A Mg 


Model space for 64 Zn 


spherical 
orbit 


sp 
energy 


spherical 
orbit 


sp 
energy 




-4.03 


lPl/2 


-8.3 


04/2 


0.08 


1P3/2 


-10.2 


04/2 


-5.00 


0/5/2 


-8.8 






0/ 7 /2 


-14.4 






0#9/2 


-4.4 
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Figure 1: Variation of energy with temperature and the quantum number 
M. The points A and B on the horizontal dashed line corresponding to E = 
-47.5 MeV indicate the change in temperature with M changing from 2 to 3. 

Figure 2: Temperature versus the inverse level density parameter K using 
eq. (18) for J = 0, 2, 4 and 6. The level density parameer a decreases with 
increase in temperature as well as angular momentum. 

Figure 3: The J^-cranking results for level density as a function of energy 
for a few values of angular momenta, J = 0, 2, 4 and 6. Fluctuations in Inp 
at J = 6 occur due to the numerical difficulty for adjusing u in r.h.s. of eq. 
(13) up to a reasonable accuracy. For example, the points a and b correspond 
to the angular momentum J a = 6.15 and J& = 5.85, respectively. 

Figure 4: Level density of 24 Mg as a function of energy with its angular 
momentum dependence extracted using eq. (14). 

Figure 5: Same as figure 4 with orientation fluctuation correction. 

Figure 6: Inverse level density parameter K as a function of temperature at 
fixed values of angular momenta, J = 0, 8, 16 and 28 for 64 Zn. 

Figure 7: Level density of 64 Zn as a function of energy at fixed values of 
angular momenta obtained in J^-cranking. 

Figure 8: Same as figure 7, with angular momentum dependence extracted 
using eq. (14). 

Figure 9: Level density as a function of excitation energy E* (obtained from 
eqs. (17) and (18)) for fixed values of angular momenta. 



19 



